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Abstract 

Tunable oscillatory modes of electric-field domains in doped semiconductor 
superlattices are reported. The experimental investigations demonstrate the 
realization of tunable, GHz frequencies in GaAs-AlAs superlattices covering 
the temperature region from 5 to 300 K. The orgin of the tunable oscillatory 
modes is determined using an analytical and a numerical modeling of the 
dynamics of domain formation. Three different oscillatory modes are found. 
Their presence depends on the actual shape of the drift velocity curve, the 
doping density, the boundary condition, and the length of the superlattice. 
For most bias regions, the self-sustained oscillations are due to the formation, 
motion, and recycling of the domain boundary inside the superlattice. For 
some biases, the strengths of the low and high field domain change periodically 
in time with the domain boundary being pinned within a few quantum wells. 
The dependency of the frequency on the coupling leads to the prediction of a 
new type of tunable GHz oscillator based on semiconductor superlattices. 
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I. INTRODUCTION 



Since the observation and explanation of the Gunn-effect in the early 60 's it has been 
known that traveling field inhomogeneities are generated by negative differential velocity 
(NDV).0i The properties of these traveling Gunn-domains have been studied extensively 
theoretically and experimentallyi'S leading to several proposals for applications. The most 
important application is probably the frequency generator. Although oscillator devices based 
on the Gunn-effect operating in the range of 100 GHz have been reported,! Gunn-oscillators 
have not led to the expected breakthrough in technology. This is largely due to the difficulties 
in tuning the oscillation frequency. Furthermore, it is impossible to vary the drift velocity 
vs field characteristics {v{F) curve) except by using different bulk materials (e.g., GaAs, 
InGaAs or InP). 

These problems may be overcome by using an entirely different class of NDV exhibiting 
materials formed by semiconductor superlattices. In an applied electric field perpendicu- 
lar to the two-dimensional layers, several different transport mechanisms such as miniband 
transport, resonant tunneling, and real-space transfer from well to barriers [T —>■ X trans- 
port) give rise to complex v{F) curves with several regions of NDV as shown in Fig. In 
contrast to the above mentioned materials, the tunability of the barrier and well thickness 
as well as the control over the barrier height in the superlattice can be used to tailor the 
drift velocity vs field v{F) characteristics. By proper engineering, samples with several NDV 
regions can be manufactured with control over previously inaccessible features such as the 
actual shape of the NDV (cf. Fig. |l]). The frequency of the oscillations also depends on 
the number of moving charges in the superlattice, which may be controlled by varying the 
doping in the quantum wells. 

Currently, three different types of superlattice oscillations due to NDV have to be dis- 
tinguished : 

• The most prominent type of oscillations is found, when carriers within the miniband 
are accelerated beyond the Brillouin zone boundary, where their drift velocity becomes 
negative. In the absence of scattering, the electron wave packet then oscillates with 
the Bloch frequency = eFd/h, where d and h denote the superlattice period and 
Planck's constant, respectively.ol3 These so-called Bloch oscillations were predicted for 
superlattices by Esaki and Tsu in 1970.i Inspired by the prospect of THz frequency gen- 
erators based on superlattice Bloch oscillators, a long series of investigations targeting 
the miniband regime followed leading to the observation of damped Bloch oscillations 
in experiments with pulsed optical excitation a few years ago.! 

• A different type of oscillations occurs in the miniband regime, when the scattering 
times are shorter than the tunneling times. In this case transient charge accumulations 
traveling through the superlattice may lead to oscillations in the current. Experiments 
using pulsed optical excitation of carriers in superlattices with wide minibands showed 
damped oscillations with frequencies up to 20 GHz.0 More recently, reflection gain 
up to 60 GHz in doped superlattices with miniband widths of up to 80 meV was 
demonstrated.0 

• Very recently a third type of oscillations has been shown to exist in weakly coupled 
superlattices in fleld regions where tunneling into higher subbands takes place.li3tJ As 
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a consequence of domain formation, this type of oscillations may also be observed in 
field regions where the electrons tunnel between the lowest subbands. While the first 
two oscillation types have never been observed under constant bias conditions, this 
new type is shown to operate over a broad frequency range practically independent of 
the external circuit and without any external triggering. 

In this paper we investigate experimentally and theoretically the third type of oscilla- 
tions in lightly and moderately doped superlattices (SL). This type of oscillations may be 
generally observable in the NDV regime of the v{F) characteristics of such superlattices. In 
a companion paper ,0 we have shown that self-sustained oscillations may exist if the dimen- 
sionless doping parameter, u = eNDd/{eFM) (A^d is the 3D SL doping, and Fm is the field 
at the resonance we consider) lies between two values (z^i/iV, z/2) which can be calculated 
numerically from the local v{F) characteristics [N is the number of SL periods). If several 
tunneling resonances exist (cf. Fig. |l|(a) and (b)) the frequency of the oscillations increases 
with increasing index of the subband involved in the tunneling resonance. This observation 
is attributed to the larger drift velocities associated with the resonances involving higher 
subbands. Furthermore, several new oscillatory modes without the typically well-defined 
domains are shown to exist. These oscillatory modes can be tuned in frequency simply by 
changing the applied voltage. Because the oscillations are observable even at room tem- 
perature, we suggest to use superlattice oscillators based on resonant tunneling into higher 
subbands ClS db tunable source for high frequencies. 

II. EXPERIMENTAL RESULTS 

Four GaAs/AlAs superlattice structures, grown by molecular beam epitaxy, are discussed 
in this paper. The samples are hereafter referred to as 9.0/4.0 (doped and undoped), 9.0/1.5, 
and 13.3/2.7, where the first number refers to the GaAs well width and the second to the 
AlAs barrier thickness in nm (cf. Table |). The superlattices are embedded between two 
highly doped Alo.5Gao.5As contact layers with a doping density N£)=2xl0^^ cm~^ (Si for 
n~^- and Be for p+-type doping) forming n^-n-n~^ and p^-i-n^diodes. After evaporating 
AuGe/Ni contacts onto the top and substrate side and alloying them for ohmic connections, 
the samples are wet-etched into mesas of 120 /im diameter. The circular top contacts 
of 70 /im diameter leave a large part of the mesas uncovered to allow for optical access. 
All experiments are performed at 5 K in a He-flow cryostat using high-frequency coaxial 
cables with a bandwidth of 20 GHz. The time-averaged current-voltage data are recorded 
with a Keithley SMU236. The time-resolved current is detected with a Tektronix sampling 
oscilloscope CSA 803 using the GHz sampling head SD-32. 

A. Current-voltage characteristics 

Fig. shows the current- voltage (I-V) characteristics of the doped 9.0/4.0 sample in 
forward bias. Clearly visible are two plateau regions with almost constant current between 
0.5 and 4 V and 6.5 and 8.5 V. To prove that these plateaus are related to the subband 
resonances, time-of-flight (TOFl measurements were performed on an undoped 9.0/4.0 ref- 
erence sample {p'^-i-n^ diode) .113 The peak photocurrent is thereby taken as a measure for 
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the inverse transport time, which at low excitation densities is directly proportional to the 
drift velocity for homogeneous fields. Note that for the presentation in Fig. |^ the built-in 
voltage (1.5 V) of the p'^-i-n^ diode has been subtracted from the reverse bias voltage ap- 
plied to the p^-i-n^ diode and the sign has been inverted. The resonances for tunneling 
from the first into the second (Ci C2) as well as into the third (Ci C3) conduction sub- 
band are clearly observed at 6 and 15 V, respectively. The low-field transmission maximum 
(Ci ^ Ci), which in the case of a strongly coupled superlattice corresponds to miniband 
conduction, is not resolved in this sample, because the transport time is too long in this field 
range to be resolved with this method. However, it has been shown previously that even 
weakly coupled superlattices exhibit a negative differential drift velocity in this regime. § In 
addition to the Ci — > Cj resonances due to F-subbands in the GaAs wells and the exponen- 
tial nonresonant background, a pronounced step-like increase of the peak photocurrent is 
observed at about 10 V, which is assigned to enhanced transport through the lowest X-level 
in the AlAs barriers. The lowest X-states of the AlAs barriers are energetically located in 
such a way that transport through these X states becomes possible at voltages between the 
Ci —>■ C2 and Ci C3 resonances. In the drift velocity vs field characteristics, this T X 
resonance leads to a step-like increase of the drift velocity as indicated by the dotted line in 
Fig. 0(b). 

The voltages of the observed resonances can be compared to the calculated energy levels 
listed in Table |. For example, the bias voltages Ubias for F ^ F resonances are calculated 
using 

uuci - c,) = ^ , (1) 

where Ci denotes the energy of the injecting conduction subband and Cj the energy of the 
conduction subband, in which the carriers are tunneling into. refers to the number of 
periods and e to the electron charge. In the case of F ^ X resonances, however, the distance 
between the levels is reduced to half the superlattice period. Therefore, an additional factor 
of 2 has to be included, i.e., 

UUCi ^ X,) = 2 ^ (2) 

where Xj denotes the energy of the receiving state in the barrier. The voltages for the F 
resonances Ci — > C2 and Ci C3 are calculated using Eq. (|I|) resulting in 5.4 and 14.6 V for 
the 9.0/4.0 sample, respectively. Using Eq. (||), we obtain a resonance voltage of 8.2 V for 
the Ci — > Xi resonance in this sample. All values agree well with the observed values. Note 
that all observed resonances appear at shghtly higher voltages due to field inhomogeneities 
and screening effects. 

To investigate the electric field inside the superlattice, photoluminescence (PL) measure- 
ments have been performed.llilli Due to the quantum-confined Stark effect (QCSE), the PL 
from regions with strong electric fields is red-shifted with respect to the PL from low-field 
regions. The experiments reveal that in both plateau regions of this sample the PL consists 
of two peaks with their respective intensity depending on the applied bias voltage. This 
Stark-splitting of the PL signal is taken as direct evidence for the existence of two distinct 
electric-field domains, which are formed because of NDV and current conservation. For 
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the first plateau the consequence of this conservation law is depicted by the dashed line in 
Fig. |l|(b). Clearly, the electric field of the low-field domain (F^) corresponds to the Ci — >• Ci 
resonance peak, while the field strength of the high field domain (-Fr) lies on the next rising 
branch of the v{F) characteristics slightly below the Ci — > C2 resonance.llli2l Applied to the 
second plateau, the current conservation law (dashed lines in Fig. ^ shows that the high 
field domain in this voltage region forms as a consequence of resonant transDort through the 
Xi level in the barriers (Ci Xi transport). In contrast to previous work,cd where this new 
formation mechanism was considered for the first time, this measurements provide the first 
direct evidence relating X-levels in the barriers to domain formation. However, it should 
be mentioned that it is not very surprising to observe domain formation due to transport 
through the Xi level. For domain formation to occur, only a minimum of the v{F) curve is 
required in conjunction with a sufficiently high carrier density.El In that sense the subject 
of domain formation and conditions for oscillations can be discussed simply by looking at 
the v{F) curve as discussed in a companion paper. El 

In the 13.3/2.7 sample the calculated Ci Xi resonance (13.7 V) is located between 
the Ci C3 and Ci C4 resonances, which are calculated to be at 9.5 and 17.9 V, 
respectively. Consequently, the I-V characteristic exhibits three plateaus as can be seen 
in Fig. ^ The thin dashed line has been added to indicate the positions of the transport 
maxima for homogeneous fields. Unlike in the 9.0/4.0 sample, the PL in the plateau regions 
shows no Stark splitting indicating that the field distribution inside the SL is less sharply 
defined. A possible explanation of this behavior is the existence of rather broad resonances 
in this weakly doped superlattice. Howeverjthe I-V characteristics shows some clear signs of 
domain formation (e.g., current branchesi^^S) , but the differences between the electric field 
values of these domains seem to be rather small. Since the QCSE is quadratic in the electric 
field, the peaks of the two domains cannot be resolved in the PL experiments. Because 
stable electric field domains have been investigated intensively over the last years, we will 
not discuss the stationary regime any further. We will rather focus on the non- stationary 
behavior instead. 



B. Oscillations in different plateaus 

The absence of a regular pattern of current branches in the plateaus is a first hint that 
the field distribution inside the SL may not always be stable. Indeed, the current, e.g., in the 
lower plateau of the 9.0/4.0 sample (Fig. shows a strongly time-dependent component, 
which has been attributed to an oscillatory instability of the field distribution inside the 
superlattice.00 During one cycle of the oscillation a monopole forms near the cathode, 
drifts through the superlattice towards the anode and dissolves, while a new monopole is 
formed.^! This type of oscillation is denoted the first oscillation mode. 

It is interesting to note that the displayed trace also contains a number of small current 
spikes. Supported by simulations, voltage turn-on measurements have shown that these 
small spikes are due to well-to-well hopping of the domain boundary moving through the 
superlattice.ii Thus, the number of current spikes indicates how many superlattice periods 
are involved in the oscillation. For the displayed oscillation this leads to a number close 
to thirty, which allows to speculate that for this particular voltage the domain boundary 
crosses about three quarters of the superlattice during one cycle of the oscillation. 
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To analyze the oscillations in more detail, we have measured the voltage dependence of 
the oscillation frequency within the plateaus. The Fourier transformlll exhibits a series of 
frequency maxima, of which we take the fundamental frequency as the oscillation frequency 
of the field distribution. 



1. Voltage independent frequency 

The voltage dependence of the oscillation frequency in the two plateaus of the 9.0/4.0 
sample is plotted in Fig. § showing that the observed frequencies in the first plateau are 
weakly voltage dependent, while in the second plateau there is no dependence on the applied 
voltage at all. The noticeable increase of the frequency with voltage at the end of the first 
plateau is due to the charge monopole having less room to travel from the beginning of the 
superlattice to the position of the high field domain. In both plateaus there are regions, 
where sinusoidal oscillations could not be observed. In these regions the current contains 
a time dependent component, which is not periodic and therefore difficult to measure. In 
particular the discontinuities in the frequency at the beginning of the first plateau suggest 
that the oscillations may become chaotic, when the charge monopole oscillates near the con- 
tacts. While chaotic behavior has been shown to exist theoretically in the case of an external 
driving frequency, the boundary regions of the plateaus are still under investigation in 
this respect. Very recently driven and undriven chaos have been observed in this system.ii 
Note that chaotic behavior has also been observed at the onset of trap-dominated Gunn os- 
cillations in ultra-pure p-Ge.izl This behavior is due to intermittent switching between small 
and large amplitude current oscillations,^ although its theoretical interpretation remains 
uncertain. ii 



2. Voltage dependent frequency 

In contrast to the 9.0/4.0 sample, the frequencies of the oscillations in the 13.3/2.7 sample 
exhibit a pronounced voltage dependence as shown in Fig. |^ for three different temperatures 
in the third plateau. The oscillations begin to appear close to the resonance maximum and 
decrease monotonically with increasing voltage. Furthermore, in this sample oscillations are 
only detected near the beginning of the plateaus in the region of NDV. This observation not 
only applies to the third plateau, but also to the first and second plateaus. Thus, there is a 
clear difference between the oscillations observed in these two samples. Unlike in the 9.0/4.0 
sample, the time-resolved PL does not show an oscillating two-domain solutionjll but merely 
an intensity modulated broad line. These two differences are taken as a signature for the 
existence of a different oscillatory mode in this sample, which is not related to well-defined 
electric field domains. 

In Fig. 1^ the experimental results on the 9.0/1.5 sample are shown. In this case we used a 
25 fim mesa in order to keep the current at a reasonable level. Due to the stronger coupling, 
GHz oscillations are already observed in the lowest plateau. In addition to a strong voltage 
dependence of the frequency, there is also an abrupt jump to lower values. Higher plateaus 
could not be measured due to the large increase in current at about 8 V. Time-resolved 
PL spectra on a 120 /xm mesa do not show two well-separated lines. However, intensity 
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oscillations with the same frequencies are clearly observed. Current oscillations on these 
large mesas are also measured using pulsed electric fields to avoid the heating of the sample. 



C. Temperature dependence 

Fig. 1^ shows the current voltage characteristic of the 13.3/2.7 sample at room tempera- 
ture. From the fact that it still shows the plateau-like features related to resonant tunneling 
one can speculate that the oscillations may also be found at higher temperatures. In Fig. |^ 
the oscillations in the third plateau are observed up to 200 K, while in Fig. ^ tunable os- 
cillations indeed exist even at 295 K in the case of the second plateau. These results are 
very encouraging, because they prove that a device such as a frequency generator based 
on the recently found tunable oscillatory mode can in principle be manufactured. How- 
ever, in order to obtain well-defined resonances even at room temperature, one probably 
has to switch to another material system with stronger confinement. For example, superlat- 
tices based on InGaAs/InAlAs material system exhibit pronounced plateaus even at room 
temperature.ii Therefore, it is expected that with properly chosen parameters oscillators 
based on InGaAs/InAlAs super lattices can be operated at room temperature. 



III. THEORETICAL INTERPRETATION: THREE OSCILLATORY MODES 

In order to explain why the frequency depends in one case on the applied voltage, while 
it is voltage independent in the other case, the condition for the two oscillatory modes has to 
be investigated in more detail. Numerical calculations on the following discrete drift model 
have been performed for the electric field Fi and 3D electron density hi of the i*^ SL period 
as well as for the total current density J using the following equations, 

- Fi., = ^{fii- No) , (3) 
dF- 

e-^ + ev{F,)fi, = J , (4) 
at 

1 A $ 

Here i = 1, . . . ,N is the period index, Nd the 3D doping density, d the SL period, e the 
average permittivity, and $ the applied voltage between the two SL ends. Eqn. (||)-(|^) have 
to be supplemented using the appropriate initial conditions for the field profile Fi as well as 
with a boundary condition for Fq, which we take as 

Fi(t)-Fo(t) = c ^, (6) 

e 

where c > denotes a dimensionless charge accumulation within the first SL period due 
to the excess doping outside the SL {n'^-n-n'^ structure). We expect that the narrower 
the barriers are, the larger c > is, since the overlap between wavefunctions of adjacent 
quantum wells is larger for narrower barriers. To find a precise relation between the structure 
of the n"*" region before the SL and c, microscopic modeling of the full n~^-n-n~^ structure 
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(characterized in Ref. |20|) is needed. This is outside the scope of the present work. However, 
we can estimate the values of c by comparing experimental and calculated frequency versus 
voltage bias curves. Thus comparing Figs. |, || and with Fig. ^(b), we obtain that c 
should be very small (c < 0.001) for samples 9.0/4.0 and 13.3/2.7, whereas c should be from 
10 to 100 times larger for sample 9.0/1.5 (see the discussion below). 

The v{F) function is a phenomenological electron velocity, which is proportional to the 
tunneling probability. The main transport mechanism in weakly coupled superlattices is 
sequential resonant tunneling (SRT). Thus, v{F) has peaks at the applied voltages or electric 
field strengths corresponding to the alignment of the energy levels of adjacent quantum wells 
as given in Eq. (|1]), which are always separated by a region of NDV. The experimentally 
observed voltage dependence of the frequency in samples 9.0/4.0 and 13.3/2.7 (cf. Figs. 5 
and 6) can be explained in terms of the actual shape of the v{F) curve, the doping density 
in the wells and the charge accumulation within the first period. The crucial features in the 
v{F) curve are the separation between the peaks as well as their width. The position of the 
peaks can be approximated by Eq. (|ID. The width of the peaks and the values of the minima 
of v{F) depend on the particular scattering mechanisms, which are present in the samples. 
However, it is reasonable to assume that the peaks of the v{F) curve for the 13.3/2.7 sample 
are wider than those of the 9.0/4.0 sample, since the narrower the barriers, the wider the 
energy intervals for resonant tunneling.E3 Therefore, we will compare simulations of the 
same model using two different v{F) curves and different doping densities. We shall use 
curves v{F) with only one peak, corresponding to a Ci Ci resonance, with i = 1,2,.... 
Thus we are only dealing with the self-sustained oscillations occurring in the plateau after 
the chosen resonance. Adding more peaks does not appreciably change the self-sustained 
oscillations on each different plateau, although new features may appear for highly doped 
SLs. In fact, there are stationary nonuniform solution branches with electric field profiles 
presenting simultaneous coexistence of three or more domains if a v{F) curve with several 
peaks is considered.0'ii 

To compare experimental results belonging to three different samples to our data, it 
is convenient to render the Eqs. (|D-(|D dimensionless. Suppose we want to analyze time- 
periodic current oscillations at a certain plateau of the current-voltage characteristic. We 
adopt as the units of electric field and velocity the corresponding values at the peak of the 
velocity curve v{F) prior to the plateau, Fm and vm, respectively. Thus, we have to scale 
the magnitudes in the figures by a different factor for each sample and plateau. The time 
scale factors, St, for the 9.0/4.0 (second plateau, after the Ci C2 resonance), 13.3/2.7 
(second plateau, after the Ci C2 resonance), and 9.0/1.5 (first plateau, after the Ci Ci 
peak), samples are 3.57, 0.98, and 0.30 ns, respectively. The electric field scale factors, se, 
are on the same order of magnitude, 10^, 4.4 x 10^, and 2.4 x 10^ V/cm. We setli 

„ Fi n,- ^ J 

Ei = —— , rii 



Fm' ' Nd' cNdVm ' 

Vm a N Fm d 



where the dimensionless parameter v is defined as 

eNi 



.= ^. (8) 
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A typical self-sustained current oscillation (far from the voltage bias corresponding to 
the onset or the end of the instability) is caused by the generation, motion and annihilation 
of domain walls, which are charge monopoles.0 We want to estimate the dependence of the 
oscillation frequency or period Tp on the bias 0. There exist three different oscillation modes 
depending on the value of the bias: 

• First mode: Oscillation of the electric field profile around an almost uniform stationary 
state (Fig. |(a)). It occurs when the bias is just above the onset of the instability. It 
can be shown that the oscillation frequency decreases as the bias increases by numer- 
ically solving the problem of linear stability of the stationary electric field profile as 
discussed in a companion paper.El 

• Second mode: Monopole recycling (Fig. ^(b)), which typically occurs when the bias 
is larger. In this mode of oscillation, the period is the sum of the formation time of 
a monopole (the time it takes each new charge accumulation injected at the first SL 
period to form a monopole) plus the time elapsed until another charge accumulation 
is injected. New monopoles are formed when the current is close to its maximum 
value (/ ^ 1, or J ^ eNjjVM in dimensional units) and the monopole formation time 
is small compared to the period of the oscillation. The bias interval, over which the 
second oscillatory mode exists, can be shown to increase with sample doping or with 
the number of SL periods. 

• Third mode: Oscillation of the electric field profile about a nonuniform stationary state 
containing two electric field domains (Fig. ^(c)). This mode is sometimes present in 
weakly doped samples, which oscillate for larger biases. It then occurs near the end 
of the oscillatory branch. In this third mode of oscillation, a domain wall remains 
pinned at a given location and the field values at the low and high field domains 
oscillate in antiphase. The shape of the current oscillation is almost sinusoidal and 
the maximum current is clearly below 1 (where the accumulation layers were injected 
from the injecting contact for the mode discussed beforel. A similar oscillation mode 
has been found in simulations of imperfect superlattices.l23 

Next we report the results of numerical simulations of the self-sustained oscillations 
of the current in the discrete model for different values of the dimensionless parameters. 
The simulations yield self- sustained current-oscillations when the doping u is such that the 
corresponding stationary field profile is clearly inhomogeneous, e.g., as in Fig. |TD|. The 
current on the middle branch of stationary solutions is significantly larger than v{(j)). In 
these cases the bias has to be within a certain interval {(pa, 4>uj) {4>a > 1 is on the NDV 
branch of v{E), and (p^j may or may not be larger than the minimum Em) to generate 
self-sustained current oscillations. Both 0^ > 1 and are increasing functions of c. In 
Fig. |l3(a), the oscillatory behavior begins at = 0^ ~ 1.100 via a supercritical Hopf 
bifurcation. The amplitude of the oscillation increases with bias following a square-root law 
and the frequency is almost constant. At = 0^^ ~ 1.617 the branch of oscillatory solutions 
disappears via a second Hopf bifurcation. For longer SLs with N = 200 and the same doping. 
Fig. |Ty(e), or for a 50-well SL with larger doping, |Tn|(c), the end of the oscillatory branch 
is different: a limit cycle collides with the unstable fixed point from the middle branch 
of the Z-shaped current-voltage characteristics and disappears. This bifurcation scenario 
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results in a decrease of the frequency down to zero, when the colhsion takes place, while 
the amplitude is unchanged to lowest order. The bias interval where the oscillatory branch 
exists, 0^), shrinks as v decreases, and for v < (i/q ~ 0.073 for = 50, c = 10"''), 
there is no oscillatory solution. Close to the Hopf bifurcations the first and third modes of 
oscillation are found. We no longer see recycling and motion of domain walls. Instead, in 
all these examples, there is an interval of bistability between the self-sustained oscillations 
and the stable stationary lower branch of the Z-shaped characteristic when G (0^,0^^). 
The bistability may be hard to observe experimentally because the basin of attraction of 
the self-sustained current oscillations is very small and most initial field profiles evolve to 
stationary situations. A similar type of bistability can be found when (i) wiggles in the 
static current- voltage characteristic appear, and (ii) the upper branch of a wiggle is 



unstable against oscillatory behavior while the lower branch is stable. See Fig. 1 of Ref. |T5]for 
V = 0.3. Fig. |10] shows the variation in magnitude and frequency of the current oscillations 
as a function of the bias for different doping values and number of SL periods. 

To clarify the nature of the oscillations let us describe the second mode of oscillations, 
which occurs for most biases. Fig. |ll|(a) shows the field profiles at different times of a given 
period of the current oscillations for a 50-well SL. We can identify a field profile consisting 
of two well-formed field domains at time (1). We define the field of the domain to be the 
electric field Ei at the position where rij has a local minimum, i.e., the variation in the field 
is minimal. The domain wall between two domains is a charge monopole containing an 
excess of electrons (nj > 1). Its evolution towards the receiving contact (until it disappears) 
is captured by snapshots at times (2) to (4) in Fig. |ll](a). Let j{t) be the instantaneous 
position of the domain wall defined as the SL period for which Ui is a local maximum. Except 
for the time intervals, in which two domain walls are present simultaneously, the constant 
voltage condition implies that 

Ei{t) = E,it)~—iE,{t)-<P) , (9) 

where Ei and Eh are the values of the field in the low and high field domains, respectively. 
Within the low field domain there is a tiny inhomogeneity (contact layer) close to the 
injecting contact, which is due to the boundary condition.El As long as the low field domain 
is in the region of positive differential velocity {Ei < 1), the field profile is stable and the 
contact layer follows adiabatically the current. The position j(t) of the domain wall moves 
to the right with a certain velocity Vmon/t^, where Vmon is always significantly less than 1 
(the peak velocity in dimensionless units). Eq. @ implies that Ei, E^, and I must increase 
with time in order to fulfill the fixed voltage condition. When Ei surpasses 1, a small charge 
accumulation is injected at i = 1, which moves to the right and grows until it becomes a well- 
formed charge monopole. See the snapshots at times (3) and (4) in Fig. 0(a). The process is 
repeated in time after the old monopole has disappeared and the new one is well established. 
Note that two different monopoles are present simultaneously during a certain part of the 
oscillation period. For a long SL the displacement current is quite small during most of each 
oscillation period, and we have Ei ^ E^^\l) and E^ ^ E^^\l) [E'^^\I) < E^^\l) < E^^\l) 
are the three solutions of the equation v{E) = I for / < 1. Note that the solutions E^^^ and 
i?*^^^ disappear for I > 1]. See Fig. |T2|for a 200-period SL. Notice that now the domains are 



well-established at times for which the values of the current are different from those depicted 
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in Fig. |TT](a). Even for a smaller SL as the one in Fig. |T^, the displacement current is small, 
when the current is near its maximum value for each period [see the time intervals between 
(1) and (2) in Fig. [Tl|(b) and between (2) and (3) in Fig. |T2](b)]. Notice that the field value 
on the low field domain is near the peak of the v{E) curve during these time intervals. This 
will be useful in what follows. 



A. Dependence of the oscillation frequency with voltage 

It is important to understand the bias dependence of the frequency which could be used 
to tune the frequency of an eventual device in agreement with the experimental observations 
of Section II. The frequency might increase or decrease with bias depending on the sample 
parameters N and u, as shown in Figs. |10|(b,d,f). Let us start at the time t = to when 
El = 1 + cu [then v{Eq) = 1 > v{Ei), Ei ^ 1], the contact layer looses its stability, and a 
charge accumulation is injected at z = 1 (to ~ 129 in Fig. |T^). One period of duration Tp 
is completed at t = to + when the next charge accumulation is injected. At this instant, 
the position of the charge monopole is given by Eq. (||) with Ei = 1 and Eh ^ E^^\l) which 
yields 

Thus, the charge accumulation has to travel a distance M in the time Tp. During the 
monopole formation time T/, Ei ^ 1, and the mean velocity of the charge accumulation is 
i>~^. On the other hand, the same charge accumulation moves at a smaller mean velocity 
equal to Vmon/j^ (with Vmon < 1) for the rest of the period, Tp — Tf, once it has become a 
charge monopole. We therefore have 

M = lTf + ^{Tp-Tf). (11) 

This gives 

Tp = Tf( 1 . (12) 

We now estimate the monopole formation time Tf. The new charge accumulation Sn travels 
with a velocity of the order of v{Ei)/v ^ 1/z/ and is amplified in time via 

'^ = -v'(E,)in(t) (13) 



with an initial value Snito) ~ c. At time (3) of Fig. |T2|, to + Ty, 5n is so large that it 
separates two different well-formed field domains. Then we may consider that a new charge 
monopole has been born. This charge monopole travels at the velocity Vmon/i^ mentioned 
above and it sharpens as it moves; see the point marked by (4) in the Fig. 0(a). This 
stage lasts until the fields before and after the monopole reach the values E = E^^\l) 
and E = E^^\l), respectively, and we are back to situation (1) having completed one 
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period. Mathematically, this behavior can be well described by an asymptotic analysis in 
the continuum limit — 0, Nu = const)^ 

In order to estimate Tf, we note that it is determined by the condition 6n{tQ + 
Tf)/6n(tQ) = a where a ~ 1/c. Then Eq. (|l^) yields 

log(a) = -/ v'iEi{t))dt. (14) 

Now we obtain Ei{t) from Eq. (||) with ^ E^^\l), j{t) = M + Vmonit - to)/^^- Up to the 
first order in {t — to), this yields 

(15) 

Linearizing v'{Ei) ^ —\v" {l)\{Ei — 1), we obtain 



Tf = VuM 



2Mlog(a) ^^g^ 



\ NVmon\v"{l)\[Ei^Kl)-l] 



Now M decreases as increases, which can be seen from Eq. (|T0|). As Vmon < 1, Eq- (|T2|) 
indicates that there are two competing mechanisms influencing the dependence of Tp with 0. 
In general the monopole formation time Tf is negligible compared to Tp(0) in the following 
cases: (i) for long SLs (M in Eq. (|12D is then large), (ii) when the values of c are large (i.e., 
a small), (iii) for large doping u. In these cases, Tp(0) should be a decreasing function. We 
find that the frequency increases with (p for large values of the quantities A^, u, and c, while 
it decreases otherwise. These behaviors are illustrated in Figs. |TD| and 



We can now interpret the experimental observations according to our theoretical results. 

• Sample 9.0/4.0 has narrow peaks, high doping with respect to the others, and is short 
(N=40). Therefore, the second oscillation mode dominates for most biases with a short 
stage of monopole formation (|f"(l)| is quite large). During most of each oscillation 
period, the electric field on the low field domain takes on values close to that of the 
resonant peak. This makes it possible to distinguish two peaks in the PL data.0'0ii 
The frequency is practically independent on the bias.0 This sample is modeled by the 
velocity curve shown in the inset of Fig. p!^(a). The results of the numerical simulation 
for current oscillations and field profiles are also shown in Fig. |T^(a). 

• Sample 13.3/2.7 has wide peaks, low doping, and is also short (A^ = 50). We would 
expect the first mode of oscillation to dominate, which agrees with the experimental 
observation that the frequency is a decreasing function of the bias. In this mode the 
high field domain is never well established, and the low field domain is off-resonance 
during most parts of each oscillation period (cf. Fig. |T3(b)). This explains that no 



Stark-splitting of the PL signal is observed in the experiments for this sample. 

Sample 9.0/1.5 is similar to sample 13.3/2.7, except that its doping is 2.5 times larger 
and the SL is shorter (A^ = 40). The velocity curve used in the simulations for this 
sample is shown in the inset of Fig. 10(d). For this sample. Fig. ^ seems to indicate that 
the first two oscillation modes are important: for low bias the first mode dominates. 
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whereas for larger bias the second mode does. Fig. |^ also shows that the frequency 
increases with bias abruptly after reaching its minimum value. Numerical calculations 
suggest that this could correspond to having a larger excess charge in the first SL 
period than in the case of sample 9.0/4.0. This could be expected for a SL with 
narrower barriers (cf. Fig. 0). 

B. Self-sustained oscillations for SL with depletion contact layers 

Until now we have studied Eqs. (|) - @ with the boundary condition of Eq. (|^) and c 
being positive. Now we are going to describe what happens if —1 < c < 0, i.e., there are 
less electrons in the first well than the doping density {rii — 1 = c < 0; cf. Fig. |T^). 

As in the case of positive c, the oscillations are due to the generation, motion and anni- 
hilation of domain walls connecting domains (which are regions of almost uniform electric 
field approximately given by the zeros of / — v{E)). The difference is that now the electric 
field profile is a monotone decreasing function of the QW index: the high-field domain is 
close to the beginning of the SL and the low-field domain extends to the end of the SL. The 
domain walls are now charge-depletion layers, having less electrons than the doping density. 
Let us describe one period of the current oscillations for a long SL as the one in Fig. [T^ with 

= 200. We will assume that initially (point marked with 1 in Fig. |15|(a)) the field profile 
has two domains connected by a domain wall {E = E^^\lo) to the left of the domain wall and 
E = E^^\lo) to the right of the domain wall with an initial value of the current /q G {vm, !))• 
The domain wall is approximately centered atj/N = Y = {(f)-E^^\lo))/{E^^^{lQ)-E^^\lo)) 
so that the voltage bias condition of Eq. (^ holds. The domain wall then moves towards 
the end of the super lattice with a speed close to the instantaneous value of the current. The 
current decreases until a certain minimum value slightly smaller than Vm is reached. Then 
a new high-field domain is created (close to the beginning of the SL) and the current rises 
sharply as the two domain walls move towards the end of the SL. When the current is near 
its maximum value, the old domain wall disappears as the values of the field in the inter- 
mediate and rightmost domains coalesce. Then one period of the oscillation is completed. 
All these features can be understood by means of an asymptotic analysis to be reported 
elsewhere.ill Some of the basic ideas were sketched in Ref. ISTl. 



IV. THE SUPERLATTICE OSCILLATOR AS A TUNABLE HIGH-FREQUENCY 

GENERATOR 

Clearly, a superlattice oscillator based on the above described types of oscillations may 
be exploited in a device, e.g., by operating the superlattice as a frequency generator. In 
order to establish the use of the newly found oscillations, it is important to know what 
maximum frequency can be achieved. While experimental investigations in this direction 
are still missing, it is clear that the ultimate frequency limit will be connected with the 
maximum drift velocity, which can be achieved without destroying the shape of the v{F) 
curve. 

To first order we assume that the oscillation frequency is determined only by the transi- 
tion time of the carriers through the sample. This approach neglects a possible oscillation 
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of the domain boundary over a few periods only, but allows to give an estimate for the 
oscillation frequencies from the current density j, the 3D doping density Nj:,, and the length 
L = Nd of the sample using 

/ = ■ (17) 

In Fig. |T6| the measured oscillation frequencies are plotted versus the frequencies calculated 
using Eq. ([T7|) . For each plateau in the 9.0/4.0 sample only one measured frequency is shown, 
while for the 13.3/2.7 and 9.0/1.5 samples the measured minimum and maximum frequencies 
are plotted. The figure shows a strong correlation between drift velocity and oscillation 
frequency over several orders of magnitude. The highest observed miniband drift velocities 
are on the order of 10^ cm/ s which, as Sibille et al.S pointed out, for a superlattice of 500 nm 
length translates into a fundamental frequency of 200 GHz. However, the electric fields in 
the reported case cannot be increased arbitrarily without breaking up the miniband, while 
for the subband field-oscillator this limitation does not apply. Accordingly, the observed 
frequencies become larger for higher subband resonances, which for sufficiently thin barriers 
will push the drift velocities and the frequencies to new limits. Looking at the oscillation 
frequency as a function of barrier width, we note that a change of a factor of 2.67 in barrier 
width leads to a change in the oscillation frequency in the first plateau of a factor of about 
4000. Projecting this behavior for other material systems to feasible barrier widths, we 
conclude that the frequencies could reach the 100 GHz range. 

Hence, in contrast to Gunn-effect devices, the oscillation frequency of the superlattice 
oscillator is not only defined by the length of the sample, but also by the parameters of 
the heterostructure. In this respect the superlattice oscillator leaves much more room for 
optimization. Additionally, a frequency generator based on a superlattice oscillator with 
several plateaus may cover several orders of magnitude with one device. In particular the 
tunability within the plateaus at high frequencies allows to envision a tunable frequency 
generator based on resonant carrier transport covering the whole GHz regime, which to our 
knowledge no other device is capable of. 



V. SUMMARY AND CONCLUSIONS 

We have demonstrated that doped, weakly coupled semiconductor superlattices exhibit 
current self-oscillations with GHz frequencies, which can be tuned by simply changing the 
applied voltage. The oscillations, which are observed in GaAs-AlAs superlattices, are not 
limited to low temperatures, but in one case have been observed at room temperature. They 
occur due to an unstable electric-field domain formation. A theoretical modeling, both 
analytical and numerical, demonstrates that three different oscillatory modes can occur, 
which depend on the detailed shape of the drift velocity curve, the doping density, the 
boundary condition, and the length of the superlattice. In contrast to the Bloch-oscillator 
proposed by Esaki and Tsu,li the self-oscillations of the field distribution are in most cases due 
to a formation, motion, and recycling of a charge accumulation layer inside the superlattice. 
In some cases the domain boundary is localized over a few quantum wells, while the field 
strengths of the low and high field domain oscillate. We would like to stress that the observed 
oscillations are also different from the damped miniband oscillations observed by Le Person 
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et al.,113 because they occur only at field strengths, where the miniband is already destroyed 
by Wannier-Stark localization. 

Considered from a macroscopic point of view in the sense that a traveling charge accumu- 
lation layer generates the oscillation in the current, the oscillations in doped super lattices are 
similar to the moving field domains in Gunn diodes. In both cases the nonlinearity, which is 
responsible for the instability, is generated by negative differential velocity. The microscopic 
origin, however, is entirely different with important consequences for the frequency limit and 
the tunability. 

In conclusion, the subband superlattice oscillator has great potential for applications as a 
tunable GHz oscillator. In order to achieve such a device, it is necessary to fully understand 
the origin of the oscillations and determine all possible oscillation modes. Furthermore, 
other material systems such as InGaAs/InAlAs should be investigated in order to generate 
even higher GHz oscillations at room temperature. 
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FIGURES 



FIG. 1. (a) Conduction band profile of a superlattice with three subbands in an apphed electric 
field, (b) Drift velocity vs field characteristics of this superlattice consisting of the nonresonant 
background, the low-field maximum (Ci ^ Ci), and the first two resonant tunneling maxima. The 
dotted line indicates a possible shape of the v{F) curve in the case of aT ^ X resonance located 
between the Ci C2 and Ci C3 resonances. The dashed line is explained in the text. 

FIG. 2. Current-voltage characteristics of the doped 9.0/4.0 sample and peak photocurrent 
vs applied voltage obtained by time-of-flight experiments in the undoped 9.0/4 reference sample. 
The dashed line indicates the condition of current conservation for the second plateau assuming a 
low-field domain with a field strength corresponding to the Ci — > C2 resonance. 

FIG. 3. Current-voltage characteristics at 4 and 295 K for the 13.3/2.7 sample. The dashed 
line shows schematically the expected drift velocity vs field characteristic of this sample for a 
homogeneous field. 

FIG. 4. Current self-oscillations at 0.78 V in the first plateau of the doped 9.0/4.0 sample at 
4 K. 

FIG. 5. Voltage dependence of the oscillation frequency within the first (a) and second (b) 
plateau of the doped 9.0/4.0 sample at 4 K. The frequencies are indicated as open circles. For 
comparison the time-averaged I-V characteristics are included labeled by dots. 

FIG. 6. Oscillation frequency vs applied voltage in the third plateau of the 13.3/2.7 sample for 
three different temperatures. 

FIG. 7. Voltage dependence of the oscillation frequency in the 9.0/1.5 sample at 4 K. The 
frequencies arc indicated as open circles. For comparison the time-averaged I-V characteristic is 
included labeled by dots. 

FIG. 8. Oscillation frequency (open circles) and amplitude (squares) vs applied voltage in the 
13.3/2.7 sample at room temperature in the second plateau. 

FIG. 9. Current oscillations (left) and field profiles (right) for v = 0.1, c = 10~^, iV = 50 and 
different values of = 1.10 (a), 1.25 (b) and 1.61 (c). Note that near the extrema of the current 
oscillations, the field distribution is close to the one of the stationary state. 
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FIG. 10. (a) Stationary current- volt age characteristic (full line) with maximum and minimum 
of the oscillating current (dotted lines) N = 50. The oscillatory branch begins at (pa ~ 1.100 and 
ends at ip^ « 1.617; the interval of bistability begins at (pp f» 1.487. (b) Fundamental frequency of 
the current vs. the average electric field (bias) for N = 50. (c,d) Same as (a,b) but now u = 0.15; 
(pa 1.052, (pp Ri 1.556, (pu, ~ 1.797. (e,f) Same as (a,b) but N = 200; ~ 1-02, (pp 1.622, 
(p^ 1.868. In all cases c = 10~^. 

FIG. 11. (a) Electric field profiles at different times during one period of the current oscillation 
depicted in the inset, (b) Time evolution of the electric field values in the left, middle (when 
it exists) and right domain of the SL. The corresponding values of E^^\l{t)), i = 1,2,3, are 
represented with dashed lines. Parameter values are (p = 1.25, c = lO-'', u = 0.1, and N = 50. 

FIG. 12. Same as in Figure 11 for N = 200. 

FIG. 13. (a) Frequency vs. bias diagrams for different values of the dimensionless doping u and 
c = 0.001. (b) Frequency vs. bias diagrams for different values of c keeping u = 0.1 and N = 50 
fixed. 

FIG. 14. Current oscillations (left) and field profiles (right) simulating samples 9.0/4.0 (a) and 
13.3/2.7 (b). 

FIG. 15. (a) Current density versus time when c = —0.01. Electron density (b) and electric 
field (c) profiles during one period of the current oscillation. The numbers in (b) and (c) correspond 
to the times marked in (a). Parameter values are (p = 1.25, u = 0.35, and N = 200. 

FIG. 16. Calculated vs measured frequencies for the different plateaus in the doped samples. 
Ip, 2p, and 3p denote the first, second, and third plateau, respectively. For the 13.3/2.7 and 9.0/1.5 
samples the minimum and maximum frequencies for each plateau are shown. 
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TABLES 



TABLE I. Parameters of the four investigated samples. N denotes the number of periods and 

N2D the two-dimensional carrier density introduced by doping. Cj and Xi indicate the energy posi- 
tions of the lowest T-subbands in the wells and the lowest A-subbands in the barriers, respectively. 
The width of the lowest conduction miniband is given by Ai. The energies of the subbands mea- 
sured from the bottom of the GaAs conduction band have been obtained by conventional envelope 
function calculations. 
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